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ABSTRACT 

We explore the relationship between the hard X-ray and sub-mm populations using deep 
Chandra observations of a large, contiguous SCUBA survey. In agreement with other recent 
findings, we confirm that the direct overlap is small. Of the 17 sub-mm sources detected in 
this field at 850 /im, only one is coincident with a Chandra source. The resulting limits imply 
that the majority of SCUBA sources are not powered by AGN, unless the central engine is 
obscured by Compton-thick material with a low (< 1 per cent) scattered component. Further- 
more, since Chandra detects only ^ 5 per cent of SCUBA sources, the typical obscuration 
would need to be almost isotropic. The X-ray upper limits are so strong that in most cases 
we can also rule out a starburst SED at low redshift, suggesting that the majority of SCUBA 
sources lie at z > 1 even if they are purely starburst galaxies. Despite the low detection rate, 
we find evidence for strong angular clustering between the X-ray and sub-mm populations. 
The implication is that AGN and SCUBA sources trace the same large-scale structure but 
do not generally coincide. If bright sub-mm sources represent massive elliptical galaxies in 
formation, we suggest that (for a given galaxy) the major episode of star-formation must be 
distinct from the period of observable quasar activity. 

Key words: galaxies: active - quasars: general - X-rays: general - X-rays: galaxies - diffuse 
radiation - galaxies: evolution - galaxies: formation - galaxies: starburst 



1 INTRODUCTION 

Our understanding of the high-redshift Universe changed dramat- 
ically with the advent of the SCUBA array at the James Clerk 
Maxwell Telescope. A number of groups have since announced 
the results from deep submillimetre surveys, all of which come 
to broadly the same conclusion. It appears that a significant (per- 
haps dominant) fraction of the star-formation in the high-redshift 
Universe (z > 2) took place in highly dust-enshrouded galaxies 
(Small et al. 1997, Hughes et al. 1998, Barger et al. 1998, Bales 
et al. 1999). These exceptionally luminous systems are likely to be 
the analogues of the Ultra-Luminous Infrared Galaxies (ULIRGs) 
seen locally. The major difference, however, is their space density. 



At high redshift ULIRG-like systems dominate the cosmic energy 
budget, while locally they are rare and unusual events. The discov- 
ery of this population was heralded by many as the discovery of the 
major epoch of dust-enshrouded spheroid formation (Lilly et al. 
1999, Dunlop 2001, Granato et al. 2001). This interpretation is by 
no means a consensus, however. There have also been suggestions 
that much of the sub-mm emission may originate from an extended 
cirrus component at lower redshift (Rowan-Robinson 2000). 

On a similar timescale, our understanding of the link between 
AGN and massive galaxies has also taken a leap forward. Thanks 
largely to the Hubble Space Telescope, it has become clear that 
essentially every massive galaxy in the local Universe hosts a su- 
permassive black hole (Kormendy & Richstone 1995, Magorrian 



© 0000 RAS 



2 O. Almaini et al. 



et al. 1998). In particular, the tight relationship between the black 
hole mass and the spheroidal velocity dispersion suggests a possi- 
ble link between an early epoch of quasar activity and the formation 
of the spheroid (Gebhardt et al. 2000, Ferrarese & Merritt 2000). 
So do the two phases coincide, or is there a evolutionary sequence 
from one to the other, as appears to be the case for local ULIRGs 
(Sanders et al. 1988)? This is a major unanswered question. One 
possibility is that the peak AGN activity corresponds to a termina- 
tion of the major burst of star-formation activity producing the host 
spheroid (Silk & Rees 1998, Fabian 1999). It has also been pointed 
out that a black hole requires a significant length of time to grow 
to a sufficient size to power a quasar, which may naturally lead to 
a lag between the epoch of peak star formation and a subsequent 
quasar phase (Archibald et al. 2001b). 

There are other reasons for expecting a strong link between 
AGN and SCUBA sources. Models for the X-ray background re- 
quire the existence of a vast population of heavily absorbed AGN 
(Comastri et al. 1995). There are major uncertainties in the mod- 
els (in particular the redshift and luminosity distribution of the ob- 
scured population) but if this absorbed energy is re-radiated in the 
far infrared one would expect 10 — 20 per cent of the SCUBA 
sources to contain AGN (Almaini, Lawrence & Boyle 1999, Fabian 
& Iwasawa 1999, Gunn & Shanks 1999). The analogy with lo- 
cal ULIRGs is also instructive. At the luminosities of the SCUBA 
sources (> IQi^^Lq, assuming z ~ 2) approximately 50 per cent 
of local ULIRGs show clear evidence for AGN activity (Sanders & 
Mirabel 1996). The relative roles of AGN and starburst activity in 
heating the dust remain controversial however (Genzel et al. 1998; 
cf. Vignati et al. 1999). 

Results from the first SCUBA surveys have already given 
some indication of the likely fraction which contain AGN. The first 
SCUBA-selected source to be optically identified was found to be 
a gas-rich QSO (Ivison et al. 1998; Frayer et al. 1998). Identifying 
subsequent sub-mm sources has been exceptionally difficult, but 
optical spectroscopy and SED constraints suggest that a sizeable 
fraction (~ 30 per cent) of the bright sub-mm sources are hosts 
to AGN activity (e.g. Ivison et al. 2000, Cooray 1999). In con- 
trast, however, the first joint Chandra/SCUBA observations found 
a very small overlap between the X-ray and sub-mm populations 
(Fabian et al. 2000; Severgnini et al. 2000; Homschemeier et al. 
2000; cf. Bautz et al. 2000), although there are some indications 
that a higher AGN fraction is detectable in the deepest X-ray ob- 
servations (Alexander et al. 2001). In this paper we explore the X- 
ray/sub-mm connection using Chandra observations of the largest 
contiguous sub-mm survey to date. Being somewhat brighter than 
previous sub-mm surveys this should allow tight constraints on the 
likely AGN content. The contiguous nature of this survey is also an 
ideal match with the Chandra field of view, allowing us to investi- 
gate the clustering between these populations for the first time. 



2 THE OBSERVATIONAL DATA 
2.1 The SCUBA Survey 

The 8mly SCUBA survey is the largest extragalactic sub-mm sur- 
vey undertaken so far, covering 260 arcmin^ in two regions of sky 
(ELAIS N2 and the Lockman Hole East) to a typical rms noise level 
of cr = 2.5 mJy at 850 /im. Full details of the sub-mm observations 
can be found in Scott et al. (2002). Details of the subsequent mul- 
tiwavelength follow-up can be found in Fox et al. (2002), Lutz et 
al. (2001) and Ivison et al. (2002). A significant advantage over 



previous surveys is the source extraction method, based on a data 
reduction pipeline which produces maps with independent pixels 
and uncorrelated noise images (Serjeant et al. 2002). Using a max- 
imum likelihood technique, these enable the statistical significance 
of each peak in the image to be assessed. Combined with a scries 
of Monte-Carlo simulated maps one can also examine the effects of 
completeness, confusion and the fikely fraction of spurious sources. 

We adopt the 3.5a SCUBA catalogue of Scott et al. (2002) as 
the basis for this analysis. This corresponds to 17 sources within 
the Chandra region. 

2.2 The ELAIS Deep X-ray Survey 

The N2 field is a particularly well studied region of sky, one of 
two fields observed by Chandra from the European Large Area 
ISO Survey (ELAIS, Oliver et al. 2000). Chandra observed this 
field on 2000 August 2 for 75 ks using the 2 x 2 array of ACIS- 
I CCDs. Analysis of the X-ray data reveals 91 sources within the 
16.9 X 16.9 arcmin field, to a flux limit of 5 x 10"^® erg s"^cm"^ 
(0.5 - 8.0 keV), of which 55 lie within the SCUBA map. Approx- 
imately 90 per cent of these Chandra sources have an optical iden- 
tification to a limit of _R ~ 26. Our deep optical imaging was also 
used to secure the X-ray astrometry to an accuracy of ~ 0.5 arcsec 
rms. Full details of the X-ray catalogue, source coimts and hard- 
ness ratios can be found in Manners et al. (2002). The optical/IR 
identification, photometry and preliminary spectra can be found in 
Gonzalez-Solares et al. (2002) and Willott et al. (2001). This survey 
will soon be complemented by a deep XMM observation (150ks). 
Combined with the sub-arcsec positions of Chandra, this will allow 
us to push significantly deeper, in addition to providing meaningful 
X-ray spectra and light curves. 

The full-band Chandra image, smoothed with a 2-arcsec Gaus- 
sian, is shown in Fig. 1. The positions of the 17 SCUBA sources 
detected above a threshold of 3.5cr are overlaid. 



3 THE CHANDRA/SCUBA OVERLAP 

Of the 17 SCUBA sources detected above 3.5cr, only one 
(iV2_850.8) contains a Chandra source within a 5 arcsec error ra- 
dius (corresponding to a 90 per cent positional error at 850/im). 
This is shown in Fig. 2. Based on the density of Chandra sources 
(1200deg"^), we expect only 0.13 of the 17 SCUBA error circles 
to contain a Chandra source by chance. By Poisson statistics, this 
corresponds to an 11 per cent probability of at least one chance 
alignment (although this could be increased if the populations are 
correlated). 

The SCUBA source iV2_850.8 has a flux of 5.1 ± 1.4 mJy 
at 850 /im. The associated X-ray source has a flux of 4 ± 1 x 
10"^'^ ergs^^cm"^ (0.5 - 8.0 keV). It lies on the edge of the 
SCUBA error circle, and coincides with a marginally resolved com- 
pact optical counterpart at ii = 22.45. The ring galaxy 3 arc- 
sec from the X-ray position is also a plausible SCUBA identifica- 
tion (reminiscent of the system discovered by Soucail et al. 1999). 
Given this, and the ~ 10 per cent probability of chance alignment, 
we view this Chandra/SCUBA coincidence with some caution. We 
note, however, that our deep VLA imaging of this field reveals a 
counterpart which coincides very well with the Chandra position, 
with a flux density of 108 ± 40 /ily at 1.4 GHz (Ivison et al., in 
preparation). A photometric redshift based on griHK photome- 
try suggests a redshift of ^ ~ 1.1, with 90 per cent confidence 
Unfits in the range 0.85 < z < 1.15 (assuming either a heavily 
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Figure 1. 16.9 x 16.9 arcmin Chandra image (75 ks exposure) with the positions of 17 SCUBA sources (> 3.5o") overlaid from the 8mjy survey (with 
large, 10 arcsec radius error circles for illustrative purposes). The dashed line indicates the extent of the SCUBA coverage. The Chandra field is centred on 
16'^36™47.0= + 41°01'34" (J2000). 



reddened young starburst or an elliptical galaxy SED). The X-ray 
source is relatively hard, although fairly typical for a source at this 
faint flux (see Manners et al. 2002 for a comparison of hardness ra- 
tios). Assuming z = 1.1, the X-ray spectrum is consistent with an 
absorbing column of 5 ± 2 x lO^^cm"^ (for an intrinsic power law 
with a = 0.7). Optical spectroscopy combined with sub-mm/mm 
interferometry may ultimately reveal the true counterpart to this 
SCUBA source. 

For the remaining SCUBA sources we estimate Sa X-ray up- 
per limits of 4 - 7 X lO"^** ergs~^cm"^ (0.5 - 8 keV). Further 
details of the Chandra flux limit and its spatial dependence can be 
found in Manners et al. (2002). Accurate SCUBA positions lead 
to tighter constraints on the X-ray flux in some cases, particularly 
where our 1.4GHz VLA radio observations coincide with a plausi- 
ble optical/IR identification (Ivison et al. 2002). 



3.1 X-ray stacking analysis 

Although only 1/17 SCUBA sources are detected by Chandra, we 
performed a stacking analysis on the remaining 16 sources to esti- 
mate their mean X-ray flux. Extracting X-ray data using a 7-arcsec 
radius around each sub-mm source (corresponding to a 2a SCUBA 
error circle) we produced a Chandra image with an effective expo- 
sure time of 1.2Ms. This gave a weak detection of X-ray flux at the 
level l.l±0.4x 10"^^ erg s~^cm"^ (0.5-8.0 keV), corresponding 
to a mean flux of 6.9 ± 2.5 x 10"^'' erg s"^cm"^ (0.5 - 8.0 keV) 
per SCUBA source. 



4 COMPARISON WITH TEMPLATE SEDS 

As described in Section 3, only one of the 17 SCUBA sources co- 
incides with a Chandra source. The vast majority do not yield any 
X-ray emission, although a stacking analysis yields a measurable 
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Figure 2. Optical R-band image of the region surrounding the SCUBA source A'^2_850.8 (left) with a 5-arcsec radius SCUBA error circle. This image has a 
magnitude limit of R ~ 26. On the edge of this eiTor circle lies the faint Chandra source N2:19 (ellipse) which has an optical counterpart with R = 22.45. 
It is unclear whether this, the neighbouring merging system or a fainter unidentified system is the true SCUBA source, although VLA observations favour the 
Chandra position. The unsmoothed Chandra image is shown on the right with 0.5 arcsec pixels. 



X-ray flux. To convert these limits and detections into quantitative 
statements on the AGN activity present we compare these hmits 
with the observed SEDs of local AGN and starburst galaxies. Fig. 
3 displays the expected sub-mm (850 /im) to X-ray (2 keV) spec- 
tral index as a function of redshift for a number of template SEDs. 
The derivation of these SEDs is described below. 



4.1 The quasar template 

For the quasar SED, we start with the 60 to 2 keV spectral index 
for radio-quiet quasars, based on the multiwavelength observations 
of the Palomar-Green (PG) Bright Quasar Survey (Sanders et al. 
1989). Longward of 60 /im the far-infrared/submillimetre emis- 
sion is almost certainly dominated by thermal re-radiation from 
dust (Carleton et al. 1987, Hughes et al. 1993). To model this spec- 
tral shape, we use the best fitting parameters obtained by Hughes, 
Davies & Ward (in preparation), based on the largest collection 
of far-infrared and submillimetre observations of local AGN. They 
find that the 50 — 1300/im emission is consistent with thermal re- 
radiation from dust at a temperature of 35 — 40 K. This can be 
modeled by an isothermal grey-body curve: 



oc 



,3+/3 



exp{hv/kT) 



A > 60 pim 



(1) 



in which T is the temperature of the dust and /3 is the emissiv- 
ity index, derived by assuming that the grey-body is transparent to 
its own emission. Best fitting values were found to be /3 = 1.7±0.3 
and T = 37 ± 5 K. 

In the X-ray regime, for an unabsorbed AGN we adopt a spec- 
tral index with a = 1.0, consistent withX-ray studies of broad-line 
AGN (e.g. Almaini et al. 1996, Reeves et al. 1997). To model highly 
absorbed AGN, we simulate the effect of photoelectric absorption 
using the XSPEC spectral analysis package, with large neutral hy- 
drogen columns up to 10^'* atom cm~^. Much beyond this limit the 
material will become Compton thick. Predicting the exact flux at 



2 keV is somewhat meaningless given the faintness of the Chan- 
dra sources, so we fold the output spectra through the full Chandra 
response to predict the broad-band 0.5 — 8 keV flux. The model 
values of asx are then calculated by converting the full-band flux 
back to 2 keV, assuming a mean X-ray spectral index of a = 0.7 
(as assumed for any detected sources). 

Fig. 3 illustrates the resulting sub-mm to X-ray predictions 
as a function of redshift. These are calculated for a naked quasar 
(Nh = 0), AGN with large absorbing columns (iV^ = 10^^, 10^*) 
and finally a simple Compton-thick model in which only 1 per cent 
of the nuclear emission emerges through scattering. It is interesting 
to note that for a very large absorbing column (A^j^ ~ 10^'') the 
negative K-corrections in both the X-ray and sub-mm wavebands 
effectively cancel out, leaving a flux ratio essentially unchanged 
with redshift. 

4.2 The starburst template 

Starburst galaxies are also luminous sources of hard X-ray emis- 
sion, although typically 2 — 3 orders of magnitude lower than AGN. 
The primary source of hard X-ray emission is from short-lived mas- 
sive X-ray binaries and supernovae (Helfand et al. 2001, Natarajan 
& Almaini 2000). Studies of local starbursts indicate a tight rela- 
tionship between the far-infrared and X-ray luminosities (David, 
Jones & Forman 1992), suggesting that the X-ray luminosity is a 
strong tracer of the underlying star formation. To predict the ob- 
served sub-mm to X-ray ratio, we use the mean 60^m to X-ray ratio 
from a recent compilation of starburst SEDs obtained by Schmitt et 
al. (1997). We note that previous authors have used the spectrum 
of Arp220 as a template (Fabian et al. 2000, Severgnini et al. 2000) 
which leads to significantly lower X-ray predictions. However it 
should be stressed that Arp220 is anomalously weak in X-ray emis- 
sion compared to other luminous infrared galaxies (Iwasawa et al. 
2001) so we do not consider it to be representative. Nevertheless, 
we include predictions based on the the SED of this galaxy for 
comparison. 
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Figure 3. The observed sub-mm to X-ray spectral index (asx) as a function of redshift. The dotted curve shows the predictions for radio-quiet quasars with 
range of photoelectric absorbing columns, plus a simple Compton thick model in which only 1 per cent of the nuclear flux emerges. The dashed curves show 
the predicted tracks for a typical starburst galaxy, based on a the mean SED compiled by Schmitt et al. (1997), and (for comparison) an X-ray-quiet starburst 
such as Arp-22() . Further details of these models can be found in Section 4. The horizontal lines represent the spectral indices for the 1 detected SCUBA 
source (Af2_850.8) and the results of a stacking analysis for the remainder. Limits for individual SCUBA sources are denoted by the shaded region. 



In the X-ray regime we assume a typical spectral index with 
a = 0.7, consistent with the spectrum of well-studied local star- 
bursts (Moran, Lehnert & Helfand 1999). In the far-infrared/sub- 
mm we scale from the 60/im emission assuming the same grey- 
body spectral shape as the AGN. As noted by Lawrence (2001), 
beyond 60/im the far-infrared/sub-mm spectra of AGN and star- 
burst galaxies are essentially identical. The values of (3 and T that 
we adopt are also in good agreement with SCUBA observations of 
local infrared galaxies (Dunne et al. 2000). The resulting sub-mm 
to X-ray spectral indices as a function of redshift are illustrated in 
Fig. 3. 



4.3 A comparison 

The only SCUBA source detected by Chandra (A''2_850.8) yields a 
value of asx which is consistent with a highly absorbed AGN. An 
unabsorbed QSO spectrum is strongly ruled out. It is, however, sig- 
nificantly more X-ray luminous than expected for a starburst galaxy 
at 2 > 1. 



For the remaining SCUBA sources, an inspection of Fig. 3 
yields a number of interesting conclusions; 

(i) The SCUBA sources not detected by Chandra cannot be pow- 
ered by AGN unless the X-ray emission is completely absorbed by 
Compton thick material. For redshifts 2 < 3 any emerging scat- 
tered fraction must be lower than 1 per cent. 

(ii) All of the SCUBA sources not detected by Chandra are con- 
sistent with the predicted X-ray properties of starburst galaxies at 
z > 2. The limits on asx are in fact so strong that one can plau- 
sibly use their X-ray non-detection as a redshift constraint. A star- 
burst of this sub-mm flux at 2; < 1 should have been detected by 
Chandra. The observed local scatter on the X-ray to FIR properties 
of starbursts is not insignificant (David, Jones & Forman 1992), but 
the non-detection of > 90% of SCUBA sources suggests (at least 
statistically) that the majority must lie at 2 > 1. We note, how- 
ever, that a cooler cirrus contribution to the sub-mm luminosity 
could weaken these redshift constraints (Rowan-Robinson 2000), 
as would the prevalence of anomalously X-ray-weak objects like 
Arp220. 



© 0000 RAS, MNRAS 000, 000-000 



6 O. Almaini et al. 




Figure 4. Estimates of the 850/im flux for the 55 Chandra sources lying 
within the SCUBA map, displayed as function of X-ray flux. The mean flux 
is 1.25 ± 0.4mjy. For comparison, models are shown based on a quasar 
SED at z = 1.5 with a range of absorbing columns. 



Figure 5. Statistical cross-correlation of the 17 SCUBA sources and the 
complete sample of Chandra sources (see text). Within 100" we obtain 82 
pairs compared to 51 expected from a random distribution. This represents 
a 4.3cr excess. 



5 SCUBA LIMITS ON THE CHANDRA SOURCES 

In addition to obtaining X-ray limits on the SCUBA sources, one 
can reverse the argument and use this large SCUBA survey to ob- 
tain sub-mm limits on the 55 Chandra sources that he within the 
SCUBA map region. Based on typical AGN SEDs, we anticipate 
that most of the Chandra sources will be too faint for detection 
by SCUBA in this relatively shallow survey, but for completeness 
we estimate SCUBA fluxes at the position of each Chandra source. 
These were obtained by fitting, using the source extraction algo- 
rithm discussed in Scott et al. (2001). The results are displayed in 
Fig. 4. Only the Chandra source discussed in Section 3 yields a sig- 
nificant SCUBA detection (> 3.5a), although others are detected 
at the 2 — 3(7 level. The majority yield 3cr upper limits of < 8mJy. 
Based on our current data, these sub-mm properties are broadly 
consistent with predictions based on absorbed quasar SEDs (see 
Fig. 4). Redshifts and deeper sub-mm photometry will allow us to 
investigate these SEDs further. 

By co-adding the SCUBA beams at the positions of each 
Chandra source, we obtain a mean sub-mm flux of 1.25 ± 0.4mJy. 
A noise-weighted mean gives 0.89 ± 0.3mJy. To test the signifi- 
cance of these results, the same analysis was performed after shift- 
ing the Chandra sources to random positions within the SCUBA 
map. This was repeated 100 times. Reassuringly, a null result is 
obtained in most cases, but the distribution of mean values sug- 
gests a 1 — 2 per cent probability of producing a '3(t' detection 
by chance. Deeper sub-mm observations are clearly required to ob- 
tain significant detections of individual Chandra sources, but we 
note that these stacked results are consistent with the recent work 
of Barger et al. (2001b). 



6 THE CLUSTERING OF SCUBA AND CHANDRA 
SOURCES 

6.1 The cross-correlation 

Although the direct overlap between the SCUBA and Chandra 
sources is small, an apparent clustering signal is visible by eye (see 
Fig. 1). Much of this may be due to the 'hole' in the centre of the 
field, which is visible in both the Chandra and SCUBA data. 

To quantify this clustering signal we use a two-point cross- 
correlation function wsx {d) of sub-mm and X-ray sources. This 
was obtained by counting the number of X-ray sources in succes- 
sive annuli around the SCUBA sources. The resulting radial dis- 
tribution was then compared with the counts obtained by placing 
100, 000 random points over the Chandra field area (avoiding edges 
and the chip boundaries). wsx{0) can then be obtained as follows: 



wsx{Oi 



Nsx{e^)Nii 
NsR{e^)Nx 



(2) 



where Nr is the total number of random points, Nx the num- 
ber of X-ray sources and Nsx{9i) and NsR{0i) give the number 
of SCUBA/X-ray and SCUBA/random pairs respectively. 

Unfortunately the sensitivity of Chandra degrades with radius, 
particularly towards the corner of the chips (a factor of two by 10 
arcmin). This alters the flux limit, which may lead to a false gra- 
dient in the number density of Chandra sources. To mitigate these 
effects, we restrict the analysis to the central 7 arcmin radius from 
the Chandra focus. 

The resulting cross-correlation signal is displayed in Fig. 5, 
with Poisson error bars. A positive signal is observed at small an- 
gular separations. Within a radius of 100" we obtain 82 pairs com- 
pared to 51 expected from a random distribution. This represents a 
4.3a detection of clustering. 
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Figure 6. Showing the relative positions of the SCUBA sources (unfilled 
hexagons) and Chandra identifications. The dashed line indicates the extent 
of the sub-mm coverage. Unlike Fig. 1 , the size of the Chandra points are 
proportional to the log of their optical (R-band) flux. Filled points represent 
Chandra sources with resolved optical identifications, while the stars illus- 
trate Chandra sources with stellar counterparts. The largest filled points are 
mostly low redshift galaxies, while the smaller filled points and stars should 
be predominantly high-z AGN and quasars (1 < z < 4). 
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Figure 7. Statistical cross-coirelation of the 17 SCUBA sources and the 
'high-redshift' subset of Chandra sources (see text). 



6.2 The high-redshift Chandra sources 

A problem which plagues all sub-mm sui-veys is the extreme optical 
faintness of the sources (Smail et al. 2002). Typically most are in- 
visible even in the deepest optical imaging. For this reason, none of 
the SCUBA sources in our survey have a confirmed redshift. Never- 
theless, a combination of arguments lead to the conclusion that the 
majority must lie at 2: > 1.5. This is based on their 450/850/im 
ratios, deep VLA radio observations and their optical/IR properties 
(Smail et al. 2000, Dunlop 2000, Fox et al. 2002, Ivison et al. 2001). 
The Chandra sources, however, are known to have a broad redshift 
distribution (0 < z < 4), including the newly discovered popula- 
tion of weak AGN in low redshift spheroidal galaxies (Mushotzky 
et al. 2000). 

To investigate the clustering signal further we therefore repeat 
Chandra/SCUBA cross-correlation with a 'high-redshift' Chandra 
sample. So far we have only ~ 15 redshifts for Chandra sources 
in this field. The majority lie at z > 1.5, but we exclude the 2 
examples which lie at z < 1 (Gonzales-Solares et al. 2001). Fur- 
ther redshifts will be obtained, but in the meantime we can exclude 
many obvious low-z AGN from our sample. Formally, we exclude 
all Chandra sources with bright optical IDs (R < 22) which are 
also clearly resolved as galaxies. Based on the redshift distribution 
of galaxies at this magnitude (Cohen et al. 2000) and the identifi- 
cations from the first deep Chandra surveys (Hornschemeier et al. 
2001, Barger et al. 2001a, Tozzi et al. 2001) these should lie pre- 
dominantly al z < 1. In total this excludes 19 of the 91 Chandra 
sources. Repeating the cross-correlation with the remaining 'high- 
z' sample does not increase the amplitude of w{6) significantly 
however (see Fig. 7). Formally the significance actually drops to a 
3.8a excess within 100". 

We conclude that while high-z Chandra sources account for 
most of the cross-correlation signal, at least part of this effect may 
be produced by SCUBA sources correlated with structure at low 
redshift (z < 1). This can be seen visually in Fig. 6. Could some 
fraction of the SCUBA sources lie at low redshift, or is this caused 
by gravitational lensing? We investigate this effect further in a re- 
lated paper (Almaini et al., in preparation) by cross-correlating 
SCUBA positions with optical galaxy catalogues. 



6.3 The auto-correlation functions 

To investigate clustering within the two populations separately, we 
evaluate a two-point aM?o-correlation function. This has the same 
form as equation (2) but is based on a self-correlation within 
a given catalogue. The auto-correlation function for the SCUBA 
sources was originally evaluated in Scott et al. (2001) and is re- 
produced here. A comparison with the w{9) for the 'high-redshift' 
Chandra sources is shown in Fig. 8(a). The results are somewhat 
tentative, with neither population revealing a conclusive detection 
of clustering. 

The strength of the crcsi-correlation signal, however, suggests 
that a significant fraction of these populations are tracing the same 
structure (Section 5.1). We therefore assume that these are all 'mas- 
sive galaxies' of one form or another, combine the 'high-z' Chandra 
sources with the SCUBA sources, and evaluate the auto-correlation 
function for the joint sample. Not surprisingly, the resulting angular 
auto-correlation function is noisy, but we obtain a 3.8cr detection 
of clustering within 100 arcsec (Fig. 8b). The best fitting function, 
adopting the standard functional form w{6) = AO'"'^, gives an 
amplitude y4(l°) = 0.011 ± 0.004. This is strikingly similar to the 
amplitude of clustering for bright (K < 18) EROs (Daddi et al. 
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Figure 8. On the left we display the auto-correlation function for the SCUBA sources (circles) and 'high-redshift' Chandra sources (squares). On the right 
the auto-correlation for the combined sample is displayed. For comparison, the dashed lines shows the best-fitting correlation function for EROs (Daddi et al. 
2000). 



2000) . Given that EROs are now widely believed to be dominated 
by passively evolving elliptical galaxies at 2: = 1 — 2, if the strength 
of the Chandra/SCUBA clustering can be confirmed it may argue 
for a link between these apparently disparate populations. Could 
sub-mm sources, quasars and elliptical galaxies represent massive 
galaxies at different stages in their evolution? Based on the small 
area, the small number of sources and the tentative nature of this 
clustering signal we urge caution in jumping to premature conclu- 
sions. Furthermore, the Chandra/SCUBA population are likely to 
trace higher redshifts than the EROs from Daddi et al. (2000). Dif- 
ferent N{z) distributions could lead to them having different in- 
ferred degrees of bias, even if the angular clustering is identical. 
The strength of the clustering signal may also be boosted by grav- 
itational lensing (Almaini et al., in preparation). Nevertheless, we 
believe that these results provide motivation for a sub-mm/X-ray 
survey covering a significantly larger area. The aim will be to ex- 
ploit a very basic property of hierarchical theories for structure for- 
mation, namely the strong dependence of clustering on the halo 
mass (Baugh et al. 1998, Steidel et al. 1999, Magliocchetti et al. 

2001) . 



7 IMPLICATIONS FOR THE JOINT FORMATION OF 
SPHEROIDS AND AGN 

The major new result of this work is the detection of a strong cross- 
correlation between the Chandra and SCUBA sources. This imme- 
diately suggests that these phenomena are tracing the same large 
scale structure. We have also confirmed the findings of other recent 
authors, namely that the direct overlap between the bright SCUBA 
and Chandra populations is small (at most 10 per cent). At first 
sight, these results appear somewhat contradictory. If the SCUBA 



sources are indeed the progenitors of massive elliptical galaxies 
then they should eventually harbour exceptionally massive black 
holes (Kormendy & Richstone 1995, Magorrian et al. 1998). In par- 
ticular, the tight relationship between black hole mass and the mass 
of the spheroid (Gebhardt et al. 2000, Ferrarese & Merritt 2000) 
suggests a direct link between the epoch of spheroid formation (the 
SCUBA phase?) and the growth of the black hole. Our observa- 
tions, however, suggest that the two phenomena only coincide ~ 5 
per cent of the time. There are several plausible explanations for 
this low overlap: 

(i) The nucleus is not receiving fuel. 

With star-formation rates of ^ lOOOAf0yr~^ (as required to power 
the most luminous SCUBA sources) the resulting stellar winds and 
supernovae should certainly generate a plentiful supply of fuel. 
The transfer of this material to the sub-parsec scale of an accretion 
disk, however, is by no means guaranteed (Shlosman, Begelman & 
Frank 1990). It therefore seems plausible that some fraction of the 
SCUBA sources are simply not fueling the central black hole at the 
observed epoch, although this would appear inconsistent with the 
simplest explanation of the black-hole/spheroid mass relationship 
(i.e. a joint growth process). We also note however that at least 50 
per cent of local ULIRGs with equivalent luminosities show clear 
evidence for an active nucleus (Sanders & Mirabel 1996). 

(ii) The quasar terminates the star-formation 

A number of models have been proposed in which the onset of 
QSO activity terminates the formation of the host spheroid (Silk & 
Rees 1998, Fabian 1999, Granato et al. 2001). This would lead to a 
natural time-lag between the two phenomena. 

(iii) The black hole is growing. 

Archibald et al. (2001b) have proposed a new model in which 
the central quasar is alive but still growing. This is based on the 
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hypothesis that a black hole is likely to grow from a small seed 
(~ IOOM0) and hence, even accreting at the Eddington limit, it 
will require ~ 5 x 10* years to reach a sufficient size to power 
a quasar. At this stage the peak star-formation activity may have 
ended, leading to a natural lag between the SCUBA phase and 
the subsequent luminous quasar. An AGN with < 1 per cent of 
the final quasar luminosity would remain undetected by Chandra, 
corresponding to Eddington limited accretion onto a black hole of 
~ lO'^A/o for a typical SCUBA source at z = 3. 
(iv) The Compton-thick scenario. 

If the obscuring column is ^ 10^* atom cm~^ the continuum will 
be completely suppressed. Only nuclear emission scattered into our 
line of sight will be observed, and if this is less than ~ 1 per cent 
it will remain invisible even at the depth of our Chandra observa- 
tions (see Fig. 3). Very low scattered fractions are not uncommon 
in local AGN-dominated ULIRGs (Fabian et al. 1996). It should 
be noted, however, that a Compton-thick explanation cannot occur 
within the context of a 'Unified Scheme' as applied to local AGN. 
If all the SCUBA sources contain enshrouded AGN, the detection 
of only ~ 5% by Chandra would require very small opening angles, 
i.e. almost isotropic obscuration. This is consistent with models in 
which a nuclear starburst both fuels and obscures the active nucleus 
(Fabian et al. 1998). 

Clearly not all of these models are not mutually exclusive, e.g. 
the AGN may be obscured by Compton thick material while the 
black hole is growing. 



8 SUMMARY AND CONCLUSIONS 

We present deep Chandra observations of a wide-area SCUBA sub- 
mm survey. In agreement with other recent findings, we confirm 
that the direct overlap between the X-ray and bright sub-mm pop- 
ulations is small. Only 1/17 SCUBA sources are coincident with 
Chandra detections. Coadding the X-ray flux for the remaining 16 
sources we obtain a weak detection, but this is entirely consistent 
with starburst activity. 

By a detailed comparison with the SEDs of quasars and star- 
burst galaxies, we conclude that the majority of the bright SCUBA 
sources are not powered by AGN, unless the central engine is ob- 
scured by Compton-thick material with a low ( < 1 per cent) 
scattered fraction. It should be noted, however, that a Compton- 
thick explanation cannot occur within the context of the 'Unified 
Scheme' scenario as applied to local Seyferts. If most SCUBA 
sources contain enshrouded AGN, the detection of only 5% by 
Chandra would require a very small typical opening angle, i.e. al- 
most isotropic obscuration. 

The X-ray limits are so strong in most cases that we can also 
rule out a starburst SED at low redshift. The non-detection of 16/17 
SCUBA sources suggests (at least statistically) that the majority lie 
at z > 1 even if they are purely starburst galaxies. 

Reversing the perspective, only one of the 55 Chandra sources 
lying within the overlap region is significantly detected by SCUBA, 
although several yield detections at the 2 — 3a level. The remain- 
der show typical 3cr upper limits of < 8 mjy at 850 /im, broadly 
consistent with the expectations of absorbed AGN models. By co- 
adding the SCUBA beams at the positions of each Chandra source, 
we obtain a mean sub-mm flux of 1.25 ± 0.4mJy, consistent with 
Bargeretal. (2001b). 

Despite the absence of X-ray emission from 95 per cent of 
bright sub-mm sources, we find evidence for strong angular clus- 
tering between the Chandra and SCUBA populations (4.3(7 sig- 



nificance). The strength of this signal is consistent with the clus- 
tering seen among EROs (Daddi et al. 2000). The implication is 
that luminous AGN and SCUBA sources trace the same large scale 
structure, but /or a given massive galaxy the quasar phase and the 
peak episode of star formation do not coincide. This may be due to 
very heavy (Compton-thick) obscuration of the central nucleus, but 
could also reflect an evolutionary time-lag between the formation 
of the spheroid and the onset of (visible) quasar activity. 
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